FACILITY FORM 602 



ff 653 July 65 

AERODYNAMIC DESIGN AND ESTIMATED 
PERFORMANCE OF A TWO-STAGE 
CURTIS TURBINE FOR THE LIQUID OXYGEN 
TURBOPUMP OF THE M-l ENGINE 

By 

R. Beer 
Prepared for 

National Aeronautics and Space Administration 
Contract NAS 3-2555 



AEROJET-GENERAL CORPORATION 

SACRAMENTO, CALIFORNIA 

■ V 
\ 


MSA CR- 54764 
AGC 8800-12 


TECHNOLOGY REPORT 


AERODYNAMIC DESIGN AND ESTIMATED PERFORMANCE OF A 
TWO-STAGE CURTIS TURBINE FOR THE 
LIQUID OXYGEN TURBOPUMP OF THE M-l ENGINE 


Prepared For 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


19 November 1965 


CONTRACT NAS3-2555 


PREPARED BY; 

AEROJET-GENERAL CORPORATION 
LIQUID ROCKET OPERATIONS 
SACRAMENTO, CALIFORNIA 


AUTHOR ; R. Beer 


APPROVED; W. E. Campbell 

Assistant Manager 
M-l Turbopump Project 


TECHNICAL MANAGEMENT: 

NASA LEWIS RESEARCH CENTER 
CLEVELAND, OHIO 


TECHNICAL MANAGER; W. W. Wilcox 


APPROVED: W. F. Dankhoff 

M-l Project Manager 





ABSTRACT 


uj Z 3 .5 

/ 

A two-stage Curtis turbine was designed for use in the oxidizer turbopump 
of the M-l Engine. 

At its design point, the turbine produces 26,800 horsepower at a velocity 
ratio of .133 and an estimated efficiency of .53. 


Blunt edged turbine rotor airfoils are used throughout. Beside superior 
performance at subsonic Mach numbers , these airfoils ( in the form of hollow sheet 
metal blades) offer advantages in fabricability , thermal fatigue resistance, and 
weight savings as compared to airfoils with sharp leading and trailing edges. 
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I 


SUMMARY 


This report delineates the aero-thermodynamic design of a Curtis turbine 
designed for the oxidizer turbopump of the M-l Engine. At the design point, the 
turbine produces 26,800 horsepower at an estimated efficiency of 53 1o. 

The turbine design parameters are: 


Inlet total pressure (at nozzle inlet) 

psia 

200 

Inlet total temperature (at nozzle inlet) 

°R 

1190 

Outlet static pressure 

psia 

1.20 

Pressure ratio, total to static 

--- 

I.67 

Mass flow 

lb/ sec 

115 

Speed 

rpm 

3635 

Mean diameter of first stage 

in. 

33-00 

Blade- jet speed ratio 

— 

.133 


Blunt edged turbine rotor airfoils are used throughout . Besides superior 
performance at subsonic Mach numbers, these airfoils (in the form of hollow sheet 
metal blades) offer advantages in fabricability, thermal fatigue resistance, and 
weight savings as compared to airfoils with sharp leading and trailing edges a 

II o INTRODUCTION 


The pumping system of the liquid propellant M-l engine consists of two 
separate turbopumps, each having a direct-drive turbine a A gas generator, separate 
from the main engine, supplies the turbines, which are arranged in series, with the 
combustion products of liquid hydrogen and liquid oxygen . The gas is initially 
expanded in the fuel turbine and then further expanded in the oxidizer turbine. The 
exhaust from the oxidizer turbine is used in three heat exchangers; to heat hydrogen 
for the gimbal actuators, and to heat hydrogen and oxygen for tank pressurization * 

The oxidizer exhaust is then used to cool the lower section of the skirt of the 
main nozzle. Finally, the exhaust is ejected via a set of small nozzles to provide 
an approximate specific impulse of 260 lbf-sec (See Figure l) . 

Ibm 

Initially, single stage turbines were designed for the fuel and oxidizer 
turbopumps o The single stage turbine rotor for the liquid oxygen turbopump was 
fabricated from a solid forging and was used in the initial test series as a workhorse 
model. For the development engine a two-stage Curtis Turbine was specified and the 
aerodynamic design of this unit is the subject of this report. This turbine has not 
been tested full scale under hot conditions. However, actual experimental perform- 
ance data of the inlet manifold, turbine nozzle, and complete two-stage turbine has 
been obtained in cold air on a subscale model at Lewis Research Center. This 
effort will be reported separately by NASA. In addition, the fabrication methods 
used on both rotors and stators, and the design and fabrication of the unique 
integrated pump backpiate-turbine inlet manifold will be discussed in separate 
contractor reports. 
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Figure 1 . M-1 Engine Mockup 
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The nozzle, reversing vane and rotor assemblies 'were fabricated entirely from 
Inconel 718 component parts and weld assembled by the process of Electron Beam. 

Welding. Figures 2, 3 and 4 show photographs of the nozzle assembly (Figure 2), 
the reversing vane assembly (Figure 3) and the dual rotor assembly mounted to hub 
and shaft. (Figure 4) 

III. DESIGN 

A. REQUIREMENTS AND GAS PROPERTIES 

Figure 1 shows a photograph of the engine mockup indicating the ducting 
elements. At the design point, 14$ of the total flow available from the gas generator 
(to drive both turbines) by-passes the fuel turbine. The by-pass flow and the fuel 
turbine exhaust flow are carried through two 11.4-in. (inside diameter) cross-over 
ducts to the oxidizer turbine inlet, manifold. This manifold is considered to be part 
of the ducting and the turbine inlet conditions specified in this report apply to the 
annulus directly upstream, from the nozzles. It is assumed that the manifold provides 
uniform nozzle inlet conditions over the entire circumference. 


An engine system balance was performed based upon an assumed variation 
of the turbine efficiency versus the velocity ratio (U/ Co) as shown on Figure 5 and 
upon, estimated cross-over duct pressure losses. This engine system balance resulted 
in the design parameters presented as Table 1. 


The properties of the turbine gas, which is the combustion product of 
liquid hydrogen and liquid oxygen, are shown as Figure 6. The effect of the pressure 
upon Cp and Y is neglected. Also, Cp and Y are taken at the reference temperature 
of I.13O °R and .kept constant throughout the turbine . The gas properties used are as 
follows ; 


o/f 

.8 

n 

1.984 

"p 

r 

1.382 

r~_i 

Y 

= 2765 

i 

JL 

Y-i 

... 3.62 . 

R 

426 


:b. design philosophy 


BTU/lbm-°R 


lbf-ft 

lbm-°R. 


The maximum total- to- static efficiency of a turbine stage is pre- 
dominantly a function of the "j / Co ratio of the stage. At different IT/ Co ratios, 
maximum, efficiencies are obtained with different reaction distributions between 
stator and rotor. Figure 7 (l) shows the efficiency of a two stage Curtis turbine 
versus Co/ ' and reaction distribution. For U/Co - 0 = 133 or Co/U - 7 = 5, the 

(.1) Flugel, Gustav, Die Dampftu.rh.lnen , ihre Bereehnung und Konstruktion mlt 
einem Anhang uber die Gasturblnen , Johann Ambrosius Barth, Leipzig, 1931 > 
page 87 ” ” 
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igure 2. Nozzle Assembly 
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Figure 5. Estimated Turbine Efficiency vs Velocity Ratio 







TABLE 1 


M-l ENGINE 

LIQUID OXYGEN TURBOPUMP 
TURBINE DESIGN PARAMETERS 


Inlet, 

Total Pressure (at nozzle inlet) 

psia 

200 

Inlet 

Total. Temperature (at nozzle inlet) 

p. 

1190 

Oatle 

t Static Pressure 

psia 

120 

Press 

ure Ratio, Total- to- Static 

— 

I .67 

Maes 

Flow 

Ib/sec 

115 

Speed 


rpm 

3635 

Mean 

Diameter of First Stage 

in. 

33-00 

Blade 

-Jet Speed Ratio 

--- 

-133 
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REACTION DISTRIBUTION, 



Figure 7. Efficiency of Two-Stage Curtis Turbines vs Co/U and 

Reaction Distribution 
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reaction distribution 0-0-0 is recommended. This means that the entire static 
pressure drop (or static enthalpy drop) is taken in the nozzle of the first stage 
and no static enthalpy drops occur in the rotors of the first and second stages 
and reversing vanes. However , preliminary calculations indicated that this reaction 
distribution results in undesirably long blades in the second rotor. 

The reaction distribution finally selected, 0-3-7; has no enthalpy drop 
in the first stage rotor, jfo of over-all available enthalpy drop in the reversing 
vanes and 7 '$> of over-all available enthalpy drop in the second stage rotor. This 
results in desired blade heights with only a small loss in performance. 

C. AERO- THERMODYNAMIC DESIGN 

1 . Flow Quantities at Mean Diameter 


a. Loss Estimate 


It appeared that the rotor velocity coefficients, previously 
used at Aerojet-General in the design of turbines, were too conservative for the M-l 
Engine oxidizer turbine. This M-l oxidizer turbine has a nozzle height in excess of 
3-in,., a rotor blade width of approximately 1.5-in., and only subsonic velocities 
throughout its blading. It is expected to have high blading efficiencies. Because 
information concerning the losses in turbine bladings similar to the ones to be 
designed were not available, a published method of loss estimation was selected 
This method for loss estimation is summarized in Appendix B. 


Prior to the selection of this method, the velocity coefficient 
of the rotor blade shown in Figure 8 was calculated using five different methods for 
comparison purposes. The results were as follows: 


(3) 

Stenning' ' 
(It) 

Traupel * 

(5) 

Vavra w ' 


k r = V w 4 
• 94 
.90 
.86 


Aerojet-General (Figure 9) 

Ainley^^ (extrapolated) 

Traupel provides the most realistic answer and his method was 

above . 


.835 

.80 

selected as mentioned 


fpT Traupel, Walter, Therm isebe Turb omaschinen , Erster Band , Springer Verlag/ 

Berlin/ Got tingen/Heidelberg , 1958 

(3) Stenning, Alan H., Design of Turbines for High-Energy-Fuel-Low-Power-Output 
Applications , MIT Report No. 79 

(4) Traupel, Walter, op. cit. 

(5) Vavra, M. H., Analysis and Design of Modified 87-5 Turbine , AGLR No. 3; April 1962 

(6) Ainley, D. G. and Mathieson., G. C. R., A Method of Performance Estimation for 
Axial-Flow Turbines, R&M 2974, December 1951 
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Re c = 5-3. lO 5 
Re DK * 3.0. I0 5 
M* = .67 


Figure 8. Blade Configuration Used to Compare Different Methods 
for the Loss Prediction of a Rotor Blading 
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Figure 9. Rotor Velocity Loss Coefficient for Impulse Turbines 
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Performance 


bo 

Based upon the loss estimation discussed in the preceding sec- 
tion, the expansion process in the T-S Chart and the velocity triangles were calculated 
for the mean diameter « The pertinent results from these calculations are provided 
as Table 2. Figure 10 shows the process in the T-S Chart and Figure 11 shows the veloc- 
ity triangles at the mean diameter . A total- to- static efficiency of 56.7 °I<> was obtained 
from these calculations. This appears to represent the potential for this design 
after some development. Therefore, it was decided to use 53*0$ as a conservative 
estimate for the engine system balance (see Figure 5). These efficiencies are based 
upon the seal arrangement shown in Figure 12 with axial clearances, k , of not 
larger than <» 100-in. d 


2 9 Flow Quantities at Hub and Tip 


An untwisted sheet metal blade of constant cross-section was 
selected for all blade rows. As a result, the gas outlet angle of each blade row 
is nearly constant from hub to tip. The flow quantities at hub and tip were calcu- 
lated from radial equilibrium considerations, assuming constant gas outlet angles, 
and using the method given by Traupel(7)* The flow efficiencies were assumed to be 
constant from hub to tip. Results are presented in Table 3 and Figure 13* The loss 
coefficients used in the radial equilibrium calculations differ from those used in 
the performance calculations because the flow quantities at the hub and tip were 
calculated for a preliminary configuration which differed slightly from the final 
design selected. However, in view of the slight difference it was not considered 
necessary to recalculate the flow quantities at the hub and tip for the final 
configuration. 

D, BLADE DESIGN 
1. Solidity 

The loss system detailed in Appendix B requires that optimum solid- 
ities be used as shown on Figure B-2 of that appendix. Actual and optimum solidities 
for the four blade rows are as follows: 


Nozzle 


First Rotor Rev. Row 


Second Rotor 


Number of Blades 43 


98 


Actual 


1.54 1.525 


97 94 

1.51 1.46 


Optimum 


1.22-1, 4 2 1.47-1.80 


1.46-1.78 


1.44-1.76 


(7) Traupel, Walter, op. cit. 
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TABLE 2 


FLOW QUANTITIES AT TEE MEAN DIAMETER 



Dim 

Nozzle 

Inlet 

Nozzle 

Outlet 

Rotor I 
Outlet 

Rev. Vane 
Outlet 

Rotor II 
Outlet 

P T 

psia 

200 

193.0 

151.2 

146.1 

128.2 

P TR 

psia 

- 

174.5 

163.2 

I37.5 

134.6 

T 

T 

Or 

1190 

1190 

1133.3 

1133.3 

1099.5 

t t.r 

°R 

» 

1157.5 

1157-5 

1114.7 

1114.7 


°R 

- 

1061 

1080.3 

IO85 . 5 

1080.3 

p 

s 

psia 

- 

126.8 

126.8 

124.8 

120.0 

v 

ft /sec 

^ p 

3575 

2.290 

2172 

1388 

J 

ft./ sec 


524 

524 

524 

520 

w 

ft/ sec 


3090 

2765 

1700 

1850 

a 

degrees 

Z 90 

70.00 

-63.IO 

67.OO 

-56.50 

p 

degrees 

- 

66.65 

-68.00 

60.10 

-65.5O 

\ 

inch 

33 

33.00 

33.00 

33.00 

32.80 

h d 

inch 

- 

3.20 

3.90 

4.80 

5.55 

Ab 

inch 

- 

.400 

.500 

.4oo 

- 

M 

- 

- 

.800 

.509 

.481 

.309 

m r 


- 

.693 

.6i4 

.377 

.410 

P 

Itm/ft - 3 

- 

.0407 

.0401 

.0392 

.0379 

nu 

Ru* 

OVER- ALL 

BTJ 

1'bm 

i 

i 

- - - - - 

STAGE I 
112.2 
64.2 

L 

u 

L. 

1 

%* “ 
u/c 

_■ _o_ 

178.2 

176 

56.7/0 

.133 

STAGE II 

66.0 

48.6 

BTU/lbm 
BTL 7 / Ibm 


Number 0 

f 'blades 


Nozzle 

43 

Rotor I 

98 

Rev. Vanes 

97 

Rotor II 
94 

Blalirg 

Efficiency 


.91 

.80 

.83 

.86 

Blade 

roat Area in. 


111.5 

141.5 

184.5 

227.5 

Reaction 

(H x > 


• 9 

0 

.03 

.07 
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Figure 10. Expansion Process in T-S Chart 
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Figure 12. Seal Arrangement Assumed in Performance Calculations 
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TABLE 3 


FLOW QUANTITIES AT HUB, MEAN, AND TIP 
(Configuration slightly different than that finally selected) 




r : i 




TIP | 



Nozzle 

Inlet 

Nozzle 

Outlet 

Rotor I 
Outlet 

Rev. Vane 
Outlet 

Rotor II 
Outlet 

Nozzle 

INLET 

Nozzle 

Outlet 

Rotor I 
Outlet 

Rev. Vane 
Outlet 

Rotor II 
Outlet 

Nozzle 

INLET 

Nozzle 

Outlet 

Rotor I 
Outlet 

Rev. Vane 
Outlet 

Rotor II 
Outlet 

D 

in. 

29.79 

29.29 

29.07 

28.08 

27.18 

33.00 

33.00 

32.90 

32.90 

33.00 

36.21 

36.21 

36.73 

37.32 

38.82 

r 

ft 

1.240 

1.240 

1.210 

1.170 

1.132 

1.375 

1.375 

1.370 

1.370 

1.375 

1.510 

1.510 

1.530 

1.570 

1.618 

P T 

psia 

200.0 

_ 

- 

- 

124.0 

200.0 

- 

- 

- 

127.2 

200.0 

_ 

- 

- 

130.5 

p s 

psia 

200.0 

115.5 

119.2 

115.0 

115.0 

200.0 

126.8 

126.8 

124.8 

120.0 

200.0 

136.0 

132.2 

132.2 

124.4 

t t 

°R 

1190 

1190 

1133.8 

1133.8 

1102.5 

1190 

1190 

1133.0 

1133.0 

1100.9 

1190 

1190 

1132.4 

1132.4 

1100.3 

T 

°R 

1190 

1038 

1067.8 

1072.0 

1080.8 

1190 

1061 

1080.5 

1086.4 

1083.8 

1190 

1081 

1088.6 

1097.2 

1086.7 

e 

lb/ft 3 

- 

.0380 

.0382 

.0367 

.0364 

- 

.0407 

.0401 

.0393 

.0380 

- 

.0429 

.0415 

.0412 

.0391 

V 

ft/sec 

- 

3880 

2660 

2475 

1445 

- 

3570 

2299 

2150 

1300 

- 

3310 

2080 

1868 

1162 

u 

ft/sec 

- 

472 

461 

445 

431 

- 

524 

522 

522 

524 

- 

575 

582 

598 

615 

w 

ft/sec 

- 

3440 

2980 

2070 

1850 

- 

3080 

2770 

1680 

1760 

- 

2770 

2615 

1340 

1695 

c< 

C 

0 

70.0 

-64.1 

67.0 

-59.1 

0 

70.0 

-63.1 

67.0 

-56.6 

0 

70.0 

-61.9 

67.0 

-53. 5 

& 

O 

- 

67.3 

-68.0 

62.2 

-66.0 

- 

66.7 

-68.0 

60.0 

-66.0 

- 

65.9 

-68.0 

56.1 

-66.0 

PV r 

v X 

lbm 

ft . -sec . 

- 

63.0 

51.5 

41.4 

31.0 

- 

68.2 

57.1 

45.2 

37.4 

- 

73.1 

62.2 

48.3 

43.6 

rV- 

u 

f t .^/sec . 


4525 

-2780 

2660 

-1425 

- 

4620 

-2810 

2710 

-1488 

- 

4700 

-2810 

2685 

-1512 






Nozzle 

Inlet 


Nozzle 

Outlet 

Rotor 

Outl€ 

it 


Rev . 
Out 

Vane 

let 

Ro 

0 

tor II 
utlet 



A = 2T f t ?v 

Jr h 

r dr 

X 

- 



115.2 

115.0 


113. 

8 

114.5 








Stage I 

Stage II 



mL - 
u 

r 

2ir w p 

PV r A(rV ) dr 
'x u 




12950 




7290 





Jr h 
















General: 

H* 

= 1.415(12950 + 7290) = 28600 H? 



The calculations above are based on the following blading efficiencies: 


Nozzle 

Rotor I 

Rev. Vane 

Rotor II 

.91 

.81 

.81 

.79 



























Figure 13. Velocity Triangles at Hub, Mean, and Tip 
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Blade Pro f ile s 

Figure 14 illustrates the profiles . suggested 'by Loschge(^) and 
Traupel for . impulse bladings with large turning angles. Further 9 Traupel indicates 
that the loss prediction of Appendix B is consistent with profiles of this general 
shape o 

e The velocity coefficients of two impulse blades are compared on 

Figure 15^.9), Blade B shown on Figure 15 (blunt leading edge) has a better design 
point and off-design performance than the classical Blade A. 

The blunt profile is efficient over a large incidence range which 
guarantees good off-design performance . Possibly large pressure variations in the 
manifold which result in large variations of nozzle exit and rotor inlet velocities 
do not have a significant effect upon performance . 

It appears that sharp leading edges should be avoided whenever 
possible in a turbine having a very fast start transient and operating in the combus- 
tion products of liquid hydrogen and oxygen. Under these conditions 9 extreme heat 
transfer rates to the blades have been known to cause cracking of the sharp leading 
edges. 


Considering the small change in relative inlet and outlet angles 
between the low reaction blade rows (see Figure ll ) 9 it appears feasible to use the 
same profile for both rotors and the reversing vanes. Further , it can be noted from 
Figure 13 that there is a small, variation in the gas inlet angles from the hub to 
the tip for all blade rows. Therefore 3 an untwisted profile of constant cross- 
section having the general shape shown on Figure 14 (Loschge) was selected for both 
rotors and the reversing vanes. The variation in blade outlet angles from first to 
second rotor is achieved by an appropriate decrease in stagger angles. 

Because of the blunt leading edge ^ the axial distance between the 
stator outlet and the rotor inlet is necessarily large to permit the proper velocity 
distribution to be established at the blade nose. 


The configuration of the inlet manifold was considered in select- 
ing the nozzle profile. Because of the geometry of the turbine inlet manifold (2 
inlets at 80° included angle) * the inlet velocity is not axial in a large portion 
of the circumference; in fact both positive and negative incidence angles exist. 

The large leading edge radius and the pronounced contraction of the flow at the noz 
zle intake assure good performance at various incidence angles. 


Figures 1.6 and. 17 show the coordinates for the two 
Pitch 3 cord 3 and stagger angle at the mean diameter (33 *0C -in. ) are 
3.8 for all of the blade rows. Figure 19 shows the axial, plan under 


basic profiles, 
shown on Figure 
hot conditions. 


(87 Lcsehge* Konstrukti.onen sus dem Dampf turb inenb au ? Springer Verlag/ 
Berlin/Got t; ingen/He ide.lberg , 1955 
(9) Loschge j A,, op. cito, , Figure 23 
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Figure 14. Impulse Profiles for Large Turning Angles 
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Figure 15. Comparison of the Velocity Coefficients of Blades 

with Blunt and Sharp Leading Edges 
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Figure 17. Rotor Blade Profile and Coordinates 


SMOOTH TRANSITION! BLTWEE.M D&F’NE.D POINTS 




Figure 18. Blade Layout 




Axial Plan in Hot Condition 



3 . Determination of Cascade Exit Angies 


Four different methods were used to make an estimation 

an assumed blade 

row with the following results: 

Method 

Outlet Angle 

Ainley^ 10 ^ 

67-3 

Traupel^ 11 ) 

66.3 

Zappa*- 12 ) 

65 .O 

Markov^ 

65.1 


Because the tangential component of the velocity is responsible for 
the specific work of the stage, it is preferable to select one of the more conservative 
methods of Zappa or Markov. Zappa, is the more convenient to use; therefore, it was 
selected for the current design. This method, illustrated on Figure 20, expresses 
the gas efflux angle of a blade row as a function of the ratios, throat width to pitch 
(d/s), and trailing edge thickness to pitch, (te/s). 

4. Flow Areas 

The accurate determination of the flow areas is of the utmost importance 
for obtaining the required reaction distribution in the turbine. The performance 
calculations from the performance discussion (ill, C, 1 , b) give preliminary values 
of the required free stream blade height at the exit of the blade rows. These blade 
heights were corrected using Vavra's method(^) . The flow through a turbine blade 
row is expressed by the equation: 



(.10) Ainley, D. G., op. cit. 

( II ) Traupel, Walter, op. cit. 

(.12) Zappa, 0., Plot of Gas Efflux Angle from. Turbine Blade Cascade (Private 

Communication) based upon: MCA 'IN 3802, 1956 (Dunavant, J. C. and Erwin, 

Jo R . , I nvestigation of a Related .Series of Turbine Blade Profiles in 
Cascade" ) and MCA IN 3959/ 1957 (Dunavant, J. C., Cascade Investigation of 
a Related Series of 6- Percent 'Thick Guide Vane Profiles and Design Charts ) 

( 13 ) Markov , N. M . , Calculation of the Aerodynamic Characteristics of Turbine 
Blading , Associated Technical Services Inc., Glen Ridge, New Jersey, 1958 
(1.4) Vavra, M. E., op. cit. 
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Figure 20. Gas Efflux Angle from Turbine Blade Cascade 





with s 


n 


r 


Y " 


(r - i) 


and : 


'P 


f (Pressure Ratio, Loss) 


Eifferent loss coefficients are used for calculating efficiency and flow areas, 
Lhe mean polytropic efficiency, t] . , for the calculation of the flow area is 
obtained from the blade efficiency^with the following empirical relationships 


Nozzles rj p = 1 - .50 (1 - t] n ) 

Rotors and Reversing. Vanes; 7) = 1 - .67 (l - n ) 

ip r 

Table k shows the results of these calculations and presents a 
comparison of the blade heights obtained' with those discussed under performance 
(ill, G, 1, b)„ The above equations give the blade height at the blade throat 
while those in the performances discussion give the free stream annulus heights,, 

5 . Velocity Distribution on Profiles 

The following two methods were used to estimate the velocity 
distributions on the profiles. 

a. NASA Computer Program^ 

The accuracy of the results obtained with this method depends 
to a. large extent upon the accuracy with which the computer input is prepared. It 
proved difficult to estimate the effective channel at the inlet to the blades and 
as a result, the first points on the suction side outside of the physical channel 
were neglected in determining the diffusion parameter. Figure 21 is an example of 
the preparation of the computer input for the first rotor . 

b. Vavra Potential Flow Method 

This method yields the inlet stagnation point and velocity 
distribution for any inlet angle. Its limitations are the inherent two-dimensionality 
of the field plotter and its restriction to incompressible flow. 


(15) Stinson, ¥» Eh, Turbine Computer Program, NASA-Aerojet Computer Job No, 1009 j 
Aerojet-General CorpT Memorandum dated 27 September ”1962 

(16) Vavra, M. HI, Aero-Thermodynamics and Flow in Turbomachines , John. Wiley & Sons, 
Inc 0 , New York/ Lond on , i960 
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TABLE 4 





DETERMINATION OF 

THROAT AREAS AND 

BLADE HEIGHTS 






Nozzle 

1st Rotor 

Rev. Vanes 

2nd Rotor 

n 



- 

.91 

.80 

.83 

.86 

n P 




•955 

.867 

00 
00 
— q 

• 907 

n 



- 

1.361 

1.315 

1.324 

1.335 

P T 



psia 

200 

174.5 

151.2 

137.5 

rp 

u. 



°R 

1190 

1157.5 

1133.9 

1114.7 

P / 

s' 

P T 


- 

.634 

.727 

00 

no 

\ j \ 

.873 

&/A 

d 


lbm 

. 2 

sec-m. 

1.039 

.810 

.621 

.5025 

k 




1 

1.0055 

1.0030 

I.OO 85 

A d 

!l 

A 

. 2 
m. 

111 

141.5 

184.8 

227 .O 

h d 



in. 

3.19 

3.90 

4.82 

5.55 

h* 



in. 

3.20 

3.83 

4.78 

5.46 

h * 
d 

* 


in. 

3.20 

3.90 

4.80 

5-55 


* Appendix B Method (reference performance discussion, Section III, c, 1, b) 
** Final Selection 
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Results from the velocity distribution analysis are as follows: 


The ratio of the surface velocity to the outlet velocity is plotted 
in relationship to the distance measured along the blade surfaces from an aribtrary 
point , A (see Figure 21) (A is not the stagnation point). Because the Vavra method 
yields the inlet stagnation point, it is shown for the nozzle and the first rotor. 

The remaining two rows of blades were not investigated using the potential flow method; 
therefore, the velocity at A was assumed to be equal to the inlet velocity. 

The stagnation point on the trailing edge, B, exists only in potential 
flow. Actually, the flow will separate; therefore, this stagnation point is not shown. 
Instead it is assumed that the trailing edge surface velocities at the pressure and 
suction sides are equal to the leaving velocity downstream from the blade row. 

Two criteria were used to judge the velocity distribution: 

Diffusion Parameter D 
Kofskey^"^ defines a diffusion 

parameter as follows: 

D = D + K 

s p 

where • 

Wi 

D = 1 - 4 

s 

¥ max. 
s 

¥ min . 

D = 1 - p 


A desirable value for D is .45. It is thought that the 

suction side contributes the major part of the total losses. Therefore, diffusion 

parameters larger than .45 were accepted providing D was smaller than approximately 

s 

.25. 


(it) Kofskey, M. G., Cold-Air Performance Evaluation of a Three-Stage Turbine 

having a Blade-Jet Speed Ratio of .156 Designed for a i00,003-Pound-Tnrust 
Hydrogen-Oxygen Rocket Turbopump Application , TM-X-477.> NASA Lewis Research 
Center, Cleveland, Ohio 
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Separation Parameter P 


with 


P 



Vavra^ 


defines a separation parameter as follows: 


dS / a £- (£) 8 

(w/w^) 2.2 


l/o 


and. indicates that no separation will occur as long ass 


P^.090 (Re c ^)’ 

The separation parameter indicates that larger decelerations 
of the blade surface velocity are acceptable at the inlet to the blade than close to 
’.he trailing edge. This condition seems valid because the boundary layer builds up 
gradually along the blade surface. Apparently, a triangular velocity distribution, 
which has the maximum surface velocities near the leading edge, is optimum for the 
low reaction (rotor and reversing row) blades. 

Figure 22 (a) shows the velocity distribution obtained for 
the nozzle at the design inlet angle while Figure 22 (b) compares the velocity 
distribution of the design inlet angle with off-design conditions for incidence angles 
of 1 - + 30 -degrees. Figure indicates favorable velocity distributions; the 
selected profile seems adequate. 


The velocity distribution of the first rotor blades at 
the design point conditions are shown on Figure 23 (a), which illustrates the attempted 
triangular distribution. The values of the separation parameter are acceptable (see 
Figure 2.3 (c)). The diffusion parameter on the suction side (D ) is good. (Note 
that, the first two stations, which are outside of the physically-defined channel are 
neglected in the calculation of D,.) The diffusion parameter on the pressure side 
is large; however, it is thought 'that, the pressure side contributes only a minor 
portion to the total loss. Figures 23 (b) and 2.3 (c) are a comparison between the 
velocity distribution and the separation parameter obtained for the design point 
inlet, angle with the same parameters for incidence angles of +4 degrees and -7 degrees. 
The off-design performance of the blade seems satisfactory. 


(. 1 . 8 ) lavra, M „ H., Private Communication 
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Figure 23. Velocity Distribution and Separation Parameter for 

First Rotor, 98 Blades 
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The velocity distribution for (3 3 ~ 71-degrees in 
Figure 23 (b) (Method Vavra) closely resembles the one at the hub of Figure 23 (a) 
(Method NASA). The acceptable separation parameters for the velocity distribution 
of Figure 23 (b) p 3 - 71, indicates that velocity distributions according to 
Figure 23 (a) are acceptable even if the first two points were not neglected. 

The velocity distributions for the reversing vanes 
(93. Hades) and second rotor blades (88 blades) are shown on Figures 24 and 25, 
respectively. Again, the total diffusion factor is somewhat high because of a 
large contribution from, the pressure side. The suction side is satisfactory and 
the profiles seem adequate. 


For the second rotor, the effect of the solidity upon 
the velocity distribution was investigated for 98 , 88 , and 82 blades (see Figure 26 

The originally selected solidity (88 blades) proved to 
be too low; therefore, the number of blades for the reversing row and the second 
rotor was increased to 97 and 94 blades, respectively,, 
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Figure 24. Velocity Distribution for Reversing Vanes, 93 Blades 
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Figure 26. Velocity Distribution for Second Rotor, Comparison of 82, 

88, and 98 Blades 
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APPENDICES 



APPENDIX A 


NOMENCLATURE 




a. 


a* 


A 

b 

Ab 


c 


C 

o 

c 

p 

c 

V 

a 


D 


D 

D r 


h 


i 

J 


k r 

k 

a 



k 


p 

Local Speed of Sound a = Yg RT 

'os 

*2 2 Y 

Critical Speed of Sound a = 

Area 

Axial Blade Width 

Axial Distance Between Blade Rows 

Blade Chord 

Isentropic Velocity C^ = 2g Q J Ah' 

Specific Heat at Constant Pressure 
Specific Heat at Constant Volume 
Throat Width 
Diameter 

Diffusion Parameter 
Diffusion Parameter of Suction Side 
Diffusion Parameter of Pressure Side 
Proportionality Factor in Newton Second Law 32.17 
Blade Height 

Specific Static Enthalpy 
Incidence 

Mechanical Equivalent of Heat (778.2) 

Rotor Velocity Coefficient W./W/ 

4 4 \ 

Axial Clearance 
Radial Clearance 
Relative Roughness 
Leakage Factor 
Thermal Conductivity 


ft/sec 

ft/sec 

. 2 
in 

inch 
inch 
inch 
ft/sec 
BTU/lb- °R 
BTU/lb- °R 
inch 
inch 


Ibm-ft „ 
lbf-sec 

inch 
BTU/lb 
degree 
ft lb/BTU 


inch 

inch 


BTU 

(ft-HR- °R) 
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Length of Point on Blade Surface from Point A 
on Suction Side 

Length of Point on Blade Surface from Point A 
on Pressure Side 

Specific Work in Blading 

Specific Internal Work 

Absolute Mach Number 


Molecular Weight 

Mach Number Relative to Rotor Blade 
Mass Flow 

Number of Blades in a Row 
Folytropic Exponent 
Rotational Speed 
Mass Ratio Oxidizer/Fuel 
Separation Parameter 
Static Pressure 
Total Pressure 

Total Pressure Relative to Rotor Blade 
Prandtl-Number 

P 

Dynamic Head y - 
Gas Constant 
Radius 

Reynolds Number 

Reynolds Number based on Blade Chord 
Reynolds Number based on Hydraulic Diameter 


Degree of Reaction: R 

x 


i 1 rn r 

s2 " s4 



inch 


inch 

BTO / L Bm 


Ibm/ sec 


rpm 


psia 

psia 

psia 

psi 

lbf-ft 
ibm- °R 

ft 


Blade Fitch 


inch 



t 


Maximum Blade Thickness 

inch 

t 

e 


Trailing Edge Thickness 

inch 

T 

B 


Effective Blade Temperature 

°R 



Total Temperature 

°R 

T 

s 


Static Temperature 

°R 

T 

TR 


Total Temperature Relative to Rotor Blade 

°R 

U 


Wheel Velocity 

ft/sec 

V 


Absolute Velocity 

ft/sec 

W 


Velocity Relative to Rotor 

ft/sec 

a 

(alpha) 

Angle Between Axial Direction and Absolute 
Gas Velocity 

degrees 

P 

(beta) 

Angle Between Axial Direction and Relative 
Gas Velocity 

degrees 

r 

(gamma) 

c 

y = — 

8 C 

V 

Deviation' 

— 

6 

(delta) 

degrees 

A 

(delta) 

Prefix to' Indicate Change 

— 

c 

(zeta) 

Loss Coefficient 

— - 

% 


Blading Efficiency, total to total 

--- 

n 


Blading Efficiency, total to static 

— - 



Internal Efficiency, total to total 

-- - 

* 


Internal Efficiency, total to static 

--- 

% 

(eta) 

Turbine Efficiency Based on Total to Static 
Pressure Ratio 

— 

tf 

p 


Mean Politropic Efficiency 


T) 

{ n 


Flow Efficiency in Nozzle Blading 


T) 

i r 


Flow Efficiency in Rotor Blading 

— - 

Z 

(zeta) 

Stagger Angle 

degrees 
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e 

(theta) 

Camber Angle 

degrees 

n 

(mu) 

Absolute Viscosity 

LBm / (hr, 

p 

(rho) 

Density 

lbm/ f t3 

O’ 

(sigma) 

Blade Solidity = c/s 

--- 


SUBSCRIPTS 

1 Inlet Stator 


Outlet Stator 


3 

4 
d 
h 
m 
t 
u 

X 


Inlet Rotor 
Outlet Rotor 
At Blade Throat 
At Blade Root 
At Mean Blade Height 
At Blade Tip 
Tangential Component 
Axial Component 

Subscript preceding a symbol indicates the number of the stage 0 


SUPERSCRIPT 

3 Attached to temperature or enthalpy means value for isentropic 

expansion,, 
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APPENDIX B 


ESTIMATING LOSSES IN A .TURBINE STAGE 




Traupel^^ presents a complete and consistent system for estimating losses 
in a turbine stage. Because this loss system was used for the M-l oxidizer 
turbine design, the pertinent portion of Chapter 8.4^^ is abstracted and presented 
herein. Nomenclature deviated from that used in the report proper with respect 
to station numbers, which are defined in Figure Bl, and angles, which are defined 
in Figure B2. Those parameters -used solely in this appendix are defined as 
appropriate. 

1. Blading Efficiency 


Definition*: 


V, 


2g o^ 


= n r 


V 

Ah' + — 

n 




W ' 

vv o 




Ah' 

r 


VL 


+ 




(1) 


where: 


with: 


Hn or n. 


K = 


p 


ZUS 


l - C = 1 - (C + C + C + C ) 

v p w r zus ' 


Total loss coefficient 
Profile loss coefficient 
Wall loss coefficient 
Secondary loss coefficient 

Damping wire loss coefficient (not used for M-l turbine) 


( 2 ) 


The profile loss coefficient C is obtained from: 
P 


'p 'po /v p ■"■m ' v 6 ~m ' ~F 


with*: 


'po 


f ; “P 


X • X f, + C + C 1 

” o m 

from Figure B3 


(3) 


*For the rotor, the absolute angles should be replaced by the relative angles 
and the flow conditions at the stator exit by the flow conditions at the 
rotor exit. 

(l)Traupel, Walter, Thermische Turbomaschinen, Erster Band, Springer Verlag, 
Berlin/ Goettingen") Heidelberg, 1958, Pages 269 through 298. 
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! 



■po 

JV P 1 V / 

f — n- ; ks/c 

' 1 / 

from Figure B5 

X* - 

f (V^a*) 

from Figure B^ 

Xj, - 

f (1 - f) 

from Figure B 6 

- 

f (1 - f) 

from Figure b6 

e F - 


from Figure B7 


The wall loss coefficient K, , -which is due to the friction loss on 
hub and tip annuli, is estimated as follows: 


^w ^po ^p 


s^ sintt^ 


y f h • sina^ Y 


The second term gives the loss due to gas friction in the gap between 
stator and rotor. It usually is negligible. 


For the secondary loss coefficient C r we set: 



c = 

r 

X /ro + C s 

(5) 

th:* 

^ro = 

V, sinCt 

f \ 1 ) 

from Figure b8** 


11 

ft 

X 

X 

f (c/h) 

from Figure B9 


c = 

s 

f (h g /h) 

from Figure BIO 

(not applicable for present design) 


*See footnote on Page Bl. 

♦>ln Figure B8 the band for *» ro = f ( 9) was extrapolated. The band was 
assumed to become horizontal, similar to the loss coefficients given in 
VAVRA . 
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2. Velocity Triangle Efficiency 


2 

With Tl and T) or K and K ( T1 = K ) known, the velocity triangle 
*n r n r x 

can be calculated and with it the velocity triangle efficiencies: 


T] = Ah/Ah' 


su 


Lu 


Ah’ + 


2 2 
V o - V 2 
2g 0 D 


Ah + 


2 2 
V o " V 2 
2g o 'D 


A h’ + 


2 2 
Y o - V 2 
2g 0 3 


Static to static (7) 


Total to total 


(8) 


* 


Lu 


Ah’ + V / 2gD 

o o 


Ah + 


2 2 
V o ~ V 2 
2g 0 ^ 


Ah’ + V /2g U 

o o 


Total to static (9) 


3» Stage Efficiency 


The internal or stage efficiencies are defined as. follows: 


*See footnote on Page Bl. 
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n . 

*S1 


Ah - 2 al 
A h' 


Ah 

Ah’ 


2 Al 
A h' 


= V ~ 2 C 


T) . = T) - (C T + C T + 2^ + C + Ctj) Static to static (10) 

‘si ‘su v Ln Lr R v B ' 


With: 


Ln 

C T 

Lr 

2C 

R 


K 


Leakage loss coefficient for nozzle (stator) 

Leakage loss coefficient for rotor 

Sum of disc - and shroud friction loss coefficients 

Blade windage loss coefficient 

Moisture loss coefficient (not applicable for present design) 


As in section 2 the following additional efficiencies are defined 


% 


*u“«Ln + C Lr + ^R + ^v + 


r ) 

'B' 


Ah' 


Ah’ + 


2 2 
V o " V 2 

2g 0 ^ 


total to total (ll) 


l - n! - (c Ln + e Lr + 2C P+ c v+ eg 


“R 


Ah' 


v B ' Ah' + V o 72g 0 *3 

total to static (12) 


The leakage loss coefficients K r , C for blades without shrouds are calculated 
from: 

A 

C T = y — 

Ln 


= YL. 


Ln 


— , e T K_ 

I O . _ ’ Lr 1 

, smu. . 


A t 

Lr 


1 


sin|3. 


( 13 ) 


With' 


K I = f ( A a ) 


from Figure Bll 


•*See footnote page Bl. 
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A 


L 


Q 


Leakage area (A^ n > nozzle; A^ r > rotor) 
Annulus area=TTDh 


and for blades with shrouds from: 



^Ln 

1 

M 
(— i 

At 

Ln 

K II A Lr 


hj ^ 

n 

sina^ 

Lr ^ sin [3 g 

With*: 

Z 

= Number 

of laybrinths 



e 

= arc of 

admission 



K li 

= f ( Ah'/Vj 2 /^^) 

from Figure B12 


The loss coefficients for disc - and shroud friction are expressed by equations 

— * K 

( 15 ) and (l6). 



^R 

2.54 C M ( D h/ D m) ( D h/ h ) 
e 9 >jr 

(Disc) 


^R 

c f (W 1 * < b/h > 

9^ 

( Shroud ) 

With: 

C M 

n u. p 

- f (R ec = ) 

°m Jl 

from Figure B5 


9 

V x m 



' J «2 eSl 2 P 2 



t 

= Ah'/(U 2 /2g o ^ ) 



c f 

- f (E V \ > 

from Figure B5 


(15) 

( 16 ) 


*See footnote page Bl. 



The blade windage loss coefficient £ can be calculated from equation (17) 



il 

> 

c 1 - e ( ' 3 ° h b 

6 ¥1 e 1 ^" D m 


(17) 

With; 

c = 

.Ok + .52 h/Dm 

blades 

free, upstream 


c = 

.019 + 1.1 (.125 - h/Dm) 2 

blades 

/ 

covered, upstream 


c = 

.88-13 (h/Dm) 2 

blades 

free, downstream 


c = 

.02 + 3.0 (.125 - h/Dm) 2 

blades 

covered downstream 



Z b : 

number 

of admission segments which 




result 

in the total admission £ . 




Make if possible = 1. 


Frequently it is convenient to calculate the power loss due to friction (lO or 

ventilation (N ), rather than the loss coefficient £ ° r £ . 

V R V 


n_, = k cl. p u 3 (n Id ) 3 d 2 (18) 

R M m ' h m' h ' ' 

So 

H V ’ -5" 0 f 1 'L («) 

So 

The moisture losses being negligible in the present application are not discussed 
further . 


According to Ref. 3 the loss system described above is applicable to 
subsonic turbines having blades of the general shape of Figure B1 with a solidity 


according to Figure B2. 

© 





FIGURE B1 
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